load EcoliHyperShockEpi.mat
  
i = 5;
     figure
     hold on
     [scl,~]=size(cell_length{i});
     k = 18;
     l = 100;
     avg_length = 0;
     n = 0;
     R = 0.000555556;
     deltaT = R.*(t_vec{i}(k:l) - t_vec{i}(k));
     L0 = exp(deltaT);
      for j=1:scl
          if all(~isnan(cell_length{i}(j,k:l)))
              j
             n = n + 1;
             
            normalized_length = cell_length{i}(j,k:l)/cell_length{i}(j,k);
            avg_length = normalized_length + avg_length;
          end
      end
      avg_length = avg_length/n;
      
      
      deltaL = avg_length./L0 - 1;
      n
       plot(deltaT, deltaL)
     xlabel('Time (s)')
     ylabel('Length (\mum)')
     title(experiment{i})
     fig2pretty
